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Understanding the grain size-dependent failure behavior of polycrystaUine graphene is important for its 
appUcations both structurally and functionally. Here we perform molecular dynamics simulations to study 
the failure behavior of polycrystalline graphene by varying both grain size and distribution. We show that 
polycrystalline graphene fails in a brittle mode and grain boundary junctions serve as the crack nucleation 
sites. We also show that its breaking strength and average grain size follow an inverse pseudo Hall-Petch 
relation, in agreement with experimental measurements. Further, we find that this inverse pseudo 
Hall-Petch relation can be naturally rationalized by the weakest-link model, which describes the failure 
behavior of brittle materials. Our present work reveals insights into controlling the mechanical properties of 
polycrystalline graphene and provides guidelines for the applications of polycrystalline graphene in flexible 
electronics and nano-electronic-mechanical devices. 

The outstanding mechanical properties of pristine graphene, such as an extremely high Young's modulus of 
around 1 TPa and a high breaking strength of 130 GPa\ along with its spectacular electronic and thermal 
properties^ ^ make it a promising material for a wide range of applications in flexible electronics and nano- 
electronic-mechanical systems'*"^. However, it is still a daunting challenge to fabricate large sheets of pristine 
single-crystalline graphene. So far, large-area, monolayer graphene produced by chemical vapor deposition 
(CVD) is polycrystalline in nature and thus contains internal grain boundaries (GBs)^"^^. It is well-known that 
the properties of bulk polycrystalline materials are often dictated by the size of their grains and by the type of their 
GBs. Several groups have carried out experiments and modeling on polycrystalline graphene and reported the 
correlations between the grain size and mechanical properties, including the failure strain, the breaking strength, 
and the Young's modulus However, the conclusions drawn from these studies are often inconsistent. Some 
studies^^'^^ showed that GBs severely weaken the breaking strength of polycrystalline graphene, while other 
studies, for example, Lee et al.^^ and Rasool et al.^^ reported that its strength is only slightly reduced despite 
the presence of GBs. Undoubtedly, the understanding of grain size-dependent failure behavior of polycrystalline 
graphene is not only essential to control its mechanical properties, but also conducive for its applications in 
flexible electronics and nano-electronic-mechanical devices. 

In this paper, we perform a series of MD simulations to investigate the mechanical properties of polycrystalline 
graphene under tensile loading by varying both grain size and distribution. More specifically, we would like to 
answer the following questions: (1) What are the correlations between the average grain size and mechanical 
properties of polycrystalline graphene, such as failure strain, breaking strength, and Young's modulus? (2) How 
does the distribution of grain size affect the strength? (3) What type of defect controls the breaking strength and 
the failure pattern for polycrystalline graphene? Our present investigation reveals that the failure strength of 
polycrystalline graphene follows an inverse pseudo Hall-Petch relation and GB junctions controls the breaking 
strength of polycrystalline graphene. Furthermore, a specific type of GB junction is found to be highly prone to 
crack nucleation in polycrystalline graphene, where one of its connecting GBs is perpendicular or nearly per- 
pendicular to the loading direction. Crack preferentially starts at this specific type of GB junction, and then 
initially propagates along this GB and further branches out along other connected GBs and/or across grain 
interior, causing its ultimate failure. 

Results 

Figure l(a-c) show the top views of typical microstructures of annealed polycrystalline graphene with an average 
grain size of 3, 7 and 11 nm, respectively. Atoms are colored according to their atomic stress level. Figure 1(d) 
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Figure 1 | The atomic configurations of the annealed polycrystaUine graphene. (a-c) Top view of a 50 nm X 50 nm annealed polycrystalline graphene 
sheet with average grain size of 3, 7, and 11 nm, respectively. Atoms are colored according to their atomic stresses, (d) Side view of the annealed 
polycrystalline graphene with average grain size of 11 nm. Atoms are colored according to their out-of-plane displacement, (e-f) The close-up views of 
GBs and junctions. Atoms are colored according to their atomic stresses. 



displays the side view of the annealed polycrystalline graphene with 
average grain size of 1 1 nm, in which the atoms are colored accord- 
ing to their out-of-plane displacement. These out-of-plane deflec- 
tions are known to help minimize the energy of membrane-like 
crystalline structures^^'^^. Figure l(e-f) depict the close-up views of 
GBs and junctions, which are composed of large fractions of penta- 
gons and heptagons, and also small factions of squares, octagons, and 
vacancies. Our atomistic model for generating polycrystalline gra- 
phene follows Kotakoski et al.^^ and Zhang et al.^^ in which GBs were 
generated through a random process. A quantitative measure of the 
probabilities for the presence of squares, pentagons, heptagons, and 
octagons in polycrystalline graphene was reported by Kotakoski 
et al.^^. The GB structures generated by this model exhibit a slightly 
higher degree of disorder, typical for polycrystalline graphene grown 
by CVD^^'^^, and are consistent with previous MD simulations^^'^^"^^ 
and experimental observations^^'^^. 

The tensile stress -strain curves for polycrystalline graphene with 
an average grain size of 3, 5, 7, 9, and 11 nm are plotted in Fig. 2(a). 
Initially the gradient of the stress — strain curves increases as the 
applied strain increases, in agreement with the previous MD simu- 
lation results^^'^^. This gradient increment is a characteristic of the 
entropic elastic behavior arising from erasing the wrinkle in a thin 
membrane structure^^. As the applied strain is increased further, the 
polycrystalline graphene gradually flattens under stretching, the sp^ 
C— C bonds are directly stretched and as a result, the stress becomes 
linearly dependent on the applied strain. The slope of the stress - 
strain curve in this range gives the Young's modulus. Finally the 



polycrystalline graphene fails due to crack nucleation and propaga- 
tion. The trends of the Young's modulus, the failure strain, and the 
breaking strength taken as the peak tensile stress, as a function of the 
average grain size are plotted in Fig. 2(b-d), respectively. For each 
selected value of the average grain size, we perform simulations on 5 
randomly generated samples with the same average grain size but 
different initial grain configurations, and use the average of the 5 
simulations to consider statistical fluctuations. With increasing grain 
size, the Young's modulus increases while the failure strain decreases, 
consistent with the previous MD simulations^^'^^. For the breaking 
strength, it increases with the grain size, showing an inverse pseudo 
Hall-Petch relation. 

In order to further reveal the physical origin of this correlation 
between the average grain size and mechanical properties of poly- 
crystalline graphene, we plot the breaking strength vs. GB junction 
density in Fig. 3(a). It is seen that the breaking strength decreases 
with increasing density of GB junctions following a power law with 
an exponent of —0.05. For a brittle material, it is well-known that 
defect structure and distribution govern its failure properties. As a 
result, its failure behavior can be described by the weakest-link 
model, which states that the failure strength of a brittle material 
follows a power-law relation with the number of weak links in the 
materiaP^'^^. To apply the weakest-link model, it is necessary that the 
distribution of weak-links is broad and unbounded. In the present 
context, GB junctions are the weak links. GB junctions are formed by 
N (N > 3) GBs. The inset in Fig. 3(a) displays a representative GB 
triple junction. Each grain orientation (al, a2, and a3) is chosen 
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Figure 2 | The correlations between the grain size and mechanical properties, (a) The tensile stress-strain curves for polycrystalline graphene 

with the average grain sizes of 3, 5, 7, 9, and 1 1 nm. All simulations are performed at a strain rate of 4 X 10^ s~K (b-d) The trends of the Young's modulus, 

the failure strain, and the breaking strength as a function of the grain size, respectively. 



randomly, thus the grains can take any crystalline orientation. This 
means that the misorientation between any of the two grains can take 
any angle from 0 to 60 degrees. It is known that GBs with different 
misorientations will have different atomic structures, which in turn 
give rise to different GB energies and different failure strengths. In 
addition, the angles (Pl, P2, and p3) between any of two GBs are also 
formed randomly. These imply that there are infinitely different 
configurations of GBs, which give rise to infinitely different config- 
urations of GB junctions. Furthermore, the orientation of the GBs 
with respect to the loading direction also influences the failure. Since 
the GBs in our polycrystalline graphene are formed randomly, the 
formed GBs can take any angle with respect to the loading direction. 
Also since GB junctions with different angles with respect to the 
loading direction may have different failure strengths, therefore, they 
basically are different weak-links, and the number of such config- 
urations is also infinite. Consequently, the distribution of weak-link 
defects (GB junctions) is broad and unbounded in terms of the 
number of GB, the grain orientations, the angles between two GBs, 
and the orientation of the GBs with respect to the loading direction. 
Hence it is expected that the failure of polycrystalline graphene can 
be described by the weakest-link model and its failure strength cr^ 
should follow a power-law relation with the number of GB junctions 
or the GB junctions density pcB^ that is, Gs oc Pgb^j which is exactly 
the relation that we have observed in our MD results. The density of 
GB junctions is scaled with the grain size in an inverse quadratic 
relation, which is confirmed by our plot shown in Fig. 3(b). As a 
result, the breaking strength will also follow a power law with the 
grain size with an exponent of 0.1, and this power law fits nicely with 
the present MD simulation results as shown in Fig. 2(d), supporting 
the inverse pseudo Hall-Petch relation. Hence the present work indi- 
cates that the breaking strength and the average grain size follow an 
inverse pseudo Hall-Petch relation. 

We also consider the effect of the distribution of grain sizes on the 
mechanical properties of polycrystalline graphene. Figure 4 shows 
two typical grain size distributions of polycrystalline graphene, with 
the narrowed grain size distribution being generated from the 
Voronoi construction^°"^^, and the broad grain size distribution being 
generated from a continuous nucleation and growth construction 



forming a Johnson-Mehl microstructure^^. For both distributions, 
we maintain similar number of grains and do not vary the dimen- 
sions of the polycrystalline graphene sheet. As a result, the average 
grain size <d> and the number of GB junctions for both grain size 
distributions are maintained to be similar (note that the number of 
grains generated in a continuous nucleation and growth process is 
random and only an expected number of grains is obtained and not a 
fixed number). The breaking strength, along with the number of the 
GB junctions, for the polycrystalline graphene with both grain size 
distributions are listed in Table 1. It is noted that all the values are also 
averaged from 5 randomly generated samples with different initial 
grain configurations but the same average grain size and grain size 
distribution. Apparently, the width of the grain size distribution does 
not have any significant effect on the breaking strength. This finding 
further supports our above statement that the breaking strength of 
polycrystalline graphene is dictated by the GB junctions. 

It should be noted that the breaking strength measured experi- 
mentally also showed a grain size- dependent breaking strength. For 
example, Lee et al.^^ reported that the strength of polycrystalline 
graphene with an average grain size of —2.5 |im and 200 jim by 
means of AFM nanoindentation is 91 and 118 GPa, respectively. It 
is noted that, in their experiments, atomic force microscopy (AFM) 
nanoindentation was used to obtain the load-indentation depth 
curve. Based on this curve, the Young's modulus and breaking 
strength were then extracted\ Using the present model for polycrys- 
talline graphene, we predict that the breaking strength is 68 and 
104 GPa for a grain size of 2.5 |im and 200 |im, respectively - which 
are in good agreement with the AFM nanoindentation results. 
Hence, the relation between breaking strength vs. the grain size 
revealed in the present work provides a plausible explanation for 
experimentally measured results^^'^^. 

We have shown that the trend of the breaking strength with the 
grain size in the present work is different from those reported in the 
previous MD simulations^^'^^'^^. We show below that the reasons for 
the differences are due to the atomistic models used for polycrystal- 
line graphene. It is noted that some research groups^^'^^"^^ manually 
constructed polycrystalline graphene with regular hexagons for 
grains and only 5-7 pairs for GBs; while some other groups^^'^^ con- 
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Figure 3 | The correlations among the grain size, the density of GB 
junctions, and the breaking strength of polycrystalline graphene. (a) The 

trend of the breaking strength as a function of the GB junction density 
(Number of GB junctions per nm^). Inset is the schematic of one GB triple 
junction configuration. The grain orientations (al, a2, and a3) are 
randomly picked. Thus the misorientations between any of two grains are 
random and can take any value. In addition, the angles (pi, P2, and p3) 
between two GBs are also formed randomly. Furthermore, the triple 
junction can form any angle with respect to the loading direction. This also 
adds in more configurations for GB junctions. Overall, the number of 
different GB junctions (weak-links) is infinite, (b) The density of GB 
junctions as a function of the average grain size <d>. 

structed polycrystalline graphene with irregular polygons for grains 
and predominant 5-7 pairs together with a small fraction of squares 
and octagons for GBs. For example, in Grantab et al.'s simulations^^, 
their model structure is a bicrystalline graphene with specific 
arrangement of pentagons and heptagons along straight GBs. 
Besides, no GB junctions were included in their study. But in reality, 
there must be GB junctions in polycrystalline graphene. This issue 
was neatly addressed by Song et al.^^ In their work, grains took an 
ideal perfect hexagonal shape and GBs were manually constructed 
with ideal 5-7 pairs. A pseudo Hall-Petch relation between the failure 
strength and grain size was observed. In Kotakoski et al.'s work^^, the 
average grain size was varied from 3 to 12 nm, but a fixed number of 
grains was used for each average grain size. Their simulations showed 
that the failure strength did not depend on the average grain change 
in the range of 3-12 nm. Based on the present inverse pseudo Hall- 
Petch relation, we find that indeed the breaking strength only 
changes marginally when the average grain size varies from 3 to 
12 nm. Hence, the different results on the failure strength obtained 
by the previous studies arise from the different atomistic models used 



■ Narrowed grain size distribution 
• Broad grain size distribution 




0.5 1.0 1.5 
Grain size, dl<d> 



2.0 



Figure 4 | Two typical distributions of grain size. The polycrystaUine 
graphene with the narrowed grain size distribution is generated from the 
Voronoi construction, while the polycrystalline graphene with the broad 
grain size distribution is generated from a continuous nucleation and 
growth construction that forms a Johnson-Mehl microstructure. For these 
two grain size distributions, the number of grain and the dimensions of the 
polycrystalline graphene sheet are maintained to be similar. As a result, the 
average grain size <d> and the number of GB junctions for both grain size 
distributions are also the same. 

for polycrystalline graphene. By comparison, the atomic structure for 
polycrystalline graphene used in the present work contains randomly 
distributed grain size, shape and crystallographic orientation, thus it 
is more realistic and consistent with pervious experimental 
observations 

To further understand the failure origin and mechanism, we 
examine the deformation process of polycrystalline graphene with 
the different grain sizes under tensile loading. A typical failure pro- 
cess is shown in Fig. 5(a). We first notice that the failure process of 
polycrystalline graphene is through direct breakage of sp^ C-C bonds, 
rather than the motion of defects such as dislocations. We further 
notice that the crack preferentially starts at a GB junction (such as 
triple, quadruple or higher junctions) wherein one of the connecting 
GBs is perpendicular or nearly perpendicular to the loading dir- 
ection. We believe that this is critical to the fracture process in poly- 
crystalline graphene. After crack initiation, the crack propagates 
along this connecting GB, and then branches out either along con- 
necting GBs or across grain interior along armchair or zigzag paths. 
The close-up views of crack formation and propagation at a GB 
junction are shown in Fig. 5(b). Since GB junctions are generally 
more defected and contain a combination of 5, 7, 8-membered rings, 
they are less stiff and more compliant, and expected to have a lower 
load-carrying capacity than grain interior. As a result, the connecting 
GB which is perpendicular or near perpendicular to the loading 
direction has to carry extra load. Hence the location that connects 
this specific GB and the junction is most prone to cracking. It is noted 
that the intergranular propagation is due to the fact that GBs are 
generally weaker than the crystalline grain, and the transgranular 
cracking along armchair or zigzag directions can be explained by 
the direction- dependent edge energy, as reported in recent MD 
simulations^^ and experiments^^. 

Discussion 

We have performed MD simulations to investigate the deformation 
and failure behavior of polycrystalline graphene under uniaxial tens- 
ile loading focusing on the effect of grain size and its distribution on 
the mechanical properties. The key finding is that the breaking 
strength is controlled by GB junctions, and the failure preferentially 
starts from a specific GB junction where one of the connecting GBs is 
perpendicular or nearly perpendicular to the loading direction. We 
further show that the breaking strength and grain size follow an 
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Table 1 | The breaking strengths, along with the numbers of the GB junctions, for polycrystallinegraphene with narrow and broad grain size 
distributions 

Breaking strength (GPa) Number of GB junctions 

Narrowed grain size distribution 36.64 ± 2.004 80 

Broad grain size distribution 37.47 ± 1 .482 78 



inverse pseudo Hall -Fetch relation, which is consistent with some 
experimental measurements. We further show that the weakest-link 
model provides a physical basis for the observed inverse pseudo Hall- 
Petch relation. 

Methods 

The polycrystalline graphene samples are generated using the Voronoi tessellation 
method^°"^^. Each sample has a length of 50 nm, width of 50 nm, and average grain 
size ranging from 3 nm to 11 nm. It contains from 26 to 354 randomly oriented 
grains of various shapes and sizes. The number of grains calculated for a given volume 
depends on the choice of the average grain size. The grain positions are randomly 
distributed within the simulation sheet. To avoid grains with artificially high aspect 
ratios we ensured that the selected positions are not too close to each other. The 
created polycrystalline graphene sample is annealed to eliminate low or high-density 
regions near the GBs and junctions. We first anneal the sample at 3000 K for 50 ps 
after which the sample is quenched to 300 K during a 10 ps run allowing the sample 
to obtain its equilibrium size (pressure driven to zero)^^. 

Our MD simulations of polycrystalline graphene are performed using Large-scale 
Atomic/Molecular Massively Parallel Simulator (LAMMPS)^^. The interatomic 
interaction is described by the adaptive intermolecular reactive empirical bond order 
(AIREBO) potentiaP^. It should be noted that the smaller cutoff distance in the 
switching function of AIREBO potential should be taken in the range of 1 .92-2.0 A to 
avoid a known nonphysical post-hardening behavior^^. For the present study, this 
cutoff distance is set at 2.0 A^^. The uniaxial tensile fracture behavior of the poly- 
crystalline graphene is investigated under a constant strain rate of 4 X 10^ s"^ (slower 
than that typically used in MD simulations^^'^^'^^) in the x direction. The pressure in 
the 7 direction is adjusted to 0 GPa to satisfy the zero-stress boundary condition. A 
constant integration time step of 0.5 fs is used in all the simulations. Periodic 
boundary conditions are applied in all directions. The stress is calculated from the 
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Figure 5 | Deformation process of polycrystalline graphene. (a) Crack 
initiation and propagation of both intergranular and transgranular cracks 
in polycrystalline graphene with average grain size of 1 1 nm. The white 
dotted line indicates that the crack preferentially starts at the GB junction 
point. The red arrow indicates the transgranular crack, (b) The close-up 
views of crack formation and propagation in GB junction. Atoms are 
colored according to their atomic stress. 



normal tensor component along the loading direction of the virial stress. The atomic 
configurations are visualized with the AtomEye package^**. 
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